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We performed fluctuation analysis by means of the local scaling dimension for the strength 
function of the isoscalar (IS) and the isovector (IV) giant quadrupole resonances (GQR) in 40 Ca, 
where the strength functions are obtained by the shell model calculation within up to the 2p2h 
configurations. It is found that at small energy scale, fluctuation of the strength function almost 
obeys the Gaussian orthogonal ensemble (GOE) random matrix theory limit. On the other hand, 
\mJ ' we found a deviation from the GOE limit at the intermediate energy scale about 1.7MeV for the IS 

and at 0.9MeV for the IV. The results imply that different types of fluctuations coexist at different 
^v^j ■ energy scales. Detailed analysis strongly suggests that GOE fluctuation at small energy scale is due 

to the complicated nature of 2p2h states and that fluctuation at the intermediate energy scale is 
associated with the spreading width of the Tamm-Dancoff lplh states. 
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I. INTRODUCTION 



> . 

. Strengths of exciting high lying states in nuclei, for instance, a giant resonance, are spread over a certain energy 
interval due to the coupling to a background which is complicated and has huge degrees of freedom. As a result, the 
strength function, at a large energy scale, exhibits a global shape profile like a Lorentzian. The values of the center 
of energy or the total width depend on the properties of a probe such as multiporality, isospin, and so on, and change 
smoothly as the mass number of a nucleus. Most of these values are well understood up to now [1—3]. On the other 
i hand, at small energy scale limit comparable to the level spacing of the background states, the strength function may 
rapidly fluctuate from state to state on top of the global profile, if we neglect the escaping width due to the coupling to 
the continuum. It is believed that the fluctuation properties at this energy scale can be well simulated by the random 
matrix theory of the Gaussian orthogonal ensemble (GOE) [4-7], which has been verified by several experimental 
data: The nearest-neighbor level spacing distribution (NND) of neutron resonances near the threshold approximately 
follows the Wigner distribution [7,8]. The distribution of the reduced width (proportional to an absolute square of a 
component of the wave function) of these resonance also shows the Porter-Thomas distribution [7,8]. Furthermore, the 
NND of the Nuclear Data Ensemble shows the Wigner distribution and the A3 statistics of them shows A3 (L) cx In L 
[9]. All of these are typical signatures of the GOE. 

Thus, we know well the behavior of the strength function at both extreme small and large energy scales. From this 
knowledge, one may construct the following model to describe the strength function of highly excited states in nuclei: 
A collective state couples democratically to each of unperturbed background states, and the background Hamiltonian 
itself is the GOE random matrix. This GOE background model is closely connected with the pandemonium picture 
[10,11]. If we adjust the energy of the collective state and the coupling strength between the collective state and 
background states, the resultant strength function should be consistent with the giant resonance strength function 
at least at both extreme energy scale limits. Therefore, in order to find the difference between a realistic strength 
function and that from this model, we must investigate fluctuation properties of the strength function at intermediate 
energy scale [12,13]. 

Recently, we proposed a new method to perform fluctuation analysis of the strength function [14]. This method is 
devised to quantitatively characterize a fluctuation property as a function of energy scale by a new measure called 
the local scaling dimension (LSD), and is suitable to investigate fluctuation at intermediate energy scales. We applied 
this method to a simple doorway damping model. The doorway damping model is different from the above GOE 
background model with respect to two points: First, unperturbed background states are divided into two different 
classes so that unperturbed states in one class, which are called the doorway states, couple directly to a collective state, 
while states in another class do not. The second point is that doorway states and other unperturbed background states 
have a finite spreading width due to the mixing among them. (Note that the GOE background model corresponds to 
the large spreading width limit.) Fluctuation properties of the strength functions of both models are same at small 
energy scale limit. The result of the analysis clearly showed that fluctuation properties of both strength functions 



1 



deviate from each other at a certain intermediate energy scale, and this energy scale is closely related to the spreading 
width of the doorway states in the doorway damping model. 

In this paper, we apply the same method to a more realistic strength function of the giant resonances in a nucleus. 
Strength functions of the isoscalar (IS) and the isovector (IV) giant quadrupole resonance (GQR) in 40 Ca are obtained 
by the shell model calculation within up to the 2p2h configurations. The strength functions reproduce the experimental 
center of energy and the total width, and also show the GOE fluctuation at the small energy scale limit, as will be 
shown later. Furthermore, in this model, the doorway structure and the deviation from GOE in the background 
may be introduced in a natural way, while they are input by hand in the doorway damping model. Therefore these 
strength functions may provide the better test for the analysis. 

Note that a similar method of analysis was proposed and applied to the (e, e') or (p,p') experimental data on 208 Pb 
[15,16]. 

The paper is organized as follows: In Sec. II, we briefly explain the LSD. In Sec. Ill, we apply the method to 
the strength functions of the isoscalar (IS) and the isovector (IV) GQR. We also discuss the physical origin of the 
deviation from the GOE background model in detail. Finally, Sec. IV is devoted to conclusions. 



II. LOCAL SCALING DIMENSION 



We briefly explain the local scaling dimension in this section. See Ref. [14] for detail. 
The strength function is expressed as [17] 

S(E) = ^2S i 8(E-E i + E ), (1) 

i 

for exciting the nucleus with excitation energy E. Here Ei and Eq are the energy of discrete levels and the ground 
state energy, respectively, and Si denotes the strength of exciting the ith energy level. Let us assume strengths are 
normalized as J2i Si = 1. 

In order to quantitatively characterize how the strength function fluctuate at various energy scale, we consider 
binned distribution of the strength function S(E) by dividing whole energy interval under consideration into L bins 
with length e. Strength contained in nth bin is denoted by p n , 

Pn = Yl Si - ( 2 ) 

iGnth bin 

To characterize the distribution of the binned strengths, we introduce the moments of it, which are called in literature 
the partition function x m (e) defined by 

L 

Xm (e) = J2Pn= L (Pn)- (3) 

71=1 

Clearly, the partition functions contains, at small e limit, the information of the strength distribution , as well as the 
energy-strength correlation as a function of a bin width or an energy scale e. 

It is then helpful to employ the idea of the generalized fractal dimension [18,19], which is derived from the partition 
function and has a definite physical meaning. The concept of the generalized fractal dimension is useful only when 
a system has a multifractal nature, namely, has a self-similar structure against the change of the scale. Strength 
functions at a nucleus generally do not have such a property. Thus, we have to extend the concept of the generalized 
fractal dimension and finally reach the local scaling dimension (LSD) defined by 

D m (e) = 1 °^ Xmi£) . (4) 
m — 1 O log e 

Different from the generalized fractal dimension, the local scaling dimension characterizes "how the partition function 
scales for different energy scales" and accordingly is a function of a bin width or an energy scale e. The local 
scaling dimension, however, has an easily understandable physical meaning similar to those of the generalized fractal 
dimension: For instance, the value of D m (e) close to unity means that the strength function looks like a one dimensional 
object distributed smoothly over an energy interval when we look at the strength function at an energy scale e. On 
the other hand, it looks like a zero dimensional object, namely a dot if the value of D m (e) is close to zero. Note that 
the local scaling dimension reduces to the generalized fractal dimension when the value of D m (e) is constant over a 
whole interval of e. 
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(a) N=2048 (b) Local Scaling Dimension 




FIG. 1. (a) An example of the GOE strength function with matrix dimension of 2048. (b)The local scaling dimension D m (e) 
for the GOE strength function obtained after the ensemble average. The curves correspond to D m (e) for m = 2 to 5 from 
upper to lower. 

In actual calculation of the local scaling dimension, we define it by means of the finite difference under the change 
of a factor 2, 

D m (V2e)= 1 togX ro (2e)-log Xro (e) 
m — 1 log 2 

rather than the derivative in Eq. (4). Using the finite difference, the calculation is very simple for all the moments. 

One of the important strength functions may be the one obtained from the GOE random matrix. The local scaling 
dimension of this GOE strength function can be used as a reference to study other strength functions which are 
supposed to have a similar fluctuation properties to the GOE, as in the case of nuclear strength functions in the 
highly excited region. We show in Fig. 1 an example of the GOE strength function and the local scaling dimension 
D m (l), where I = e/d (the bin width measured in the unit of the level spacing d) is used as the scaling parameter. 

III. ANALYSIS OF THE GIANT QUADRUPOLE RESONANCES IN 40 Ca 
We apply the local scaling dimension to the strength functions of the giant quadrupole resonances in 40 Ca. 



A. Numerical Calculation of Strength Functions 

We calculated the strength functions of the isoscalar and the isovector giant quadrupole resonances in 40 Ca within 
the second Tamm-Dancoff approximation including the lplh and 2p2h excitations. Single-particle wave-functions 
and energies were obtained in terms of a Woods-Saxon potential including the Coulomb interaction. To simulate the 
bare (Hartree-Fock) single-particle energies erf, the above single-particle energies Sws were scaled by the effective 
mass to* jm as erf = Ews/( m * l m )- As the residual interaction, the Landau-Migdal-type interaction [20] including 
the density-dependence was adopted. The model space was constructed in terms of single-particle states within the 
four major shells, two below and two above the Fermi surface, and included all lplh states and 2p2h states whose 
unperturbed energies were less than 50MeV. Resultant number of lplh states and 2p2h states are 34 and 4144, 
respectively. We diagonalized the Hamiltonian within this model space and obtained the strength functions for the 
isoscalar and the isovector quadrupole operators. 

Figure 2 (a) and (b) show the strength functions for ISGQR and IVGQR, respectively. The mean energy and 
the total width (the standard deviation) are 23MeV(31MeV) and 4.0MeV(7.0MeV) for the ISGQR(IVGQR). The 
effective mass to*/to = 0.9 was used to reproduce the empirical mean energy for IVGQR in 40 Ca [3]. These strength 
functions were already used in Rcfs. [21] and [22] to study the generic properties of nuclear giant-resonance decay. 
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(a) ISGQR (b) IVGQR 
1 0.008 T 




FIG. 2. Strength functions for ISGQR (a) and IVGQR (b). 



B. Several Measures 

Before going to the detailed discussion for the local scaling dimension, we briefly show other fluctuation measures. 
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FIG. 3. The nearest neighbor level spacing distribution for levels between 20MeV and 30MeV. Level spacings were unfolded 
by the Strutinsky method with a smoothing width 5MeV. The dashed curve represents the Wigner distribution. 



Figure 3 shows the nearest- neighbor level spacing distribution (NND), where levels between 20MeV and 30MeV are 
considered (the number of levels is 804), and are unfolded by means of the Strutinsky method with a smoothing width 
5MeV [23]. The NND follows rather well the Wigner distribution, which indicates at the very small energy scale the 
level fluctuation almost obeys the GOE. We also plot in Fig. 4 the A3 statistics. We again find that at small energy 
range the A3 follows the GOE line. At L max ~ 15, however, A3 starts to deviate from the GOE line to upward. This 
value of L max roughly corresponds to the energy scale E ~ 0.18McV. 

Figure 5 shows the strength distribution where the distribution of the square-root of normalized strengths Si is 
plotted as a histogram; strengths are normalized also by Strutinsky method, detailed procedure of which will be 
explained in the next subsection. We can see that the strength distribution follows considerably well the Porter- 
Thomas distribution. 
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FIG. 4. The A3 statistics. The dashed curve represents the A3 for the GOE level fluctuation. Same as Fig. 3 for others. 
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(b) IVGQR 
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FIG. 5. The statistical distribution of square root of unfolded strengths, S i , associated with the individual levels between 
20MeV and 30MeV for (a)ISGQR and (b)IVGQR. The dashed curve represents the Porter-Thomas distribution which becomes 
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a Gaussian when plotted as a function of S i . 



C. Results of the Local Scaling Dimension 



Since we are not interested in the global shape of the strength function, it is convenient to introduce the normalized 
strength function where the global smooth energy dependence is removed from the original strength function as 
follows: 



S{Ei) 



(6) 



Here, S{E) and p(E) are obtained by averaging the calculated strength function, Eq.(l), and the level density 
p(E) = Y^i^iE — Ei + E ), respectively, by the Strutinsky method. We adopted 5MeV as a smoothing width. J\f 
is an overall normalization factor to guarantee ^ Si — 1 for the considered energy interval. Moreover, we adopt as 
energy levels the equidistant ones [14], namely, Ei — id, where d denotes average level spacing. By adopting these 
levels, we can neglect global energy dependence of the level density and the local level fluctuation. (We verified that 
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the inclusion of the local level fluctuation affects only the behavior of the LSD at small energy scales, and this effect 
can be described by the GOE.) Thus, we obtain the normalized strength function S(E); 

S{E)=H i S i S{E-E i +E ). (7) 

The local scaling dimension is derived from this normalized strength function. Finally, since the visible strengths of 
ISGQR are almost concentrated in the energy range between 20MeV and 30MeV, we adopt this energy interval for 
analysis. 



(a) ISGQR 




(b) IVGQR 
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FIG. 6. The normalized strength functions S(E) corresponding to the original ones for the ISGQR (a) and for the IVGQR 
(b) plotted in Fig. 2. A summation of strengths within an energy interval 20MeV-30MeV is normalized to unity. 



The normalized strength functions for ISGQR and IVGQR are plotted in Figs. 6 (a) and (b), respectively. 

First, let us discuss the ISGQR case in detail. The case of the IVGQR will be discussed in Sec. HID 4. Figures 
7(a) and (b) show the partition functions Xm(e) defined by Eq. (3) and the local scaling dimensions D m (e) defined by 
Eq. (5), respectively. The horizontal axes in both figures represent the bin width of energy axis in unit of d, where 
d represents the average level spacing over the energy range 20MeV and 30MeV (d ~ 12keV). From the partition 
function x m («), one sees that the fluctuation does not follow the power scaling law (or linear relation in the logdog 
plot). Correspondingly, the local scaling dimension D m (e) varies as a function of e. At small energy scale, the values 
of D m {e) gradually increase as the energy scale e increases, and almost follow the GOE line. As the energy scale 
e increases further, the values of D m (e) bend downwards and then again turn to increase. The behavior of D m (e) 
around this energy scale clearly deviates from the GOE limit. The value of the energy scale where the deviation 
is maximum is roughly ef ~ I40d ~ 1.7MeV. From these observations, we again find that at small energy scale, 
fluctuation of the strength function is essentially governed by the GOE, as one knows from other measures such as 
the NND or the strength distribution. The scaling analysis by means of the local scaling dimension, however, can 
reveal the new feature which may not be realized by other measures. Namely, the GOE fluctuation is limited to the 
small energy scale, while, at intermediate energy scales, a different fluctuation from the GOE arises. 

One may wonder that the statistical error could produce a fictitious dip, since the number of levels considered 
804 is not so large. Of course, there may be a statistical error. Accordingly, we should take the energy scale 
ef ~ I40e? ~ 1.7MeV mentioned above as a rough estimate. The existence of the dip itself, however, is not due to 
a statistical error. This would be verified from a systematic analysis in the next subsection. Note that the simple 
estimate of the standard deviation, of the mth local scaling dimension is <j(D m ) ~ y/3m/ V 'N to t [14], and in the 
present case for the second local scaling dimension, <r(D 2 ) — 0.12. 

A type of the smoothing procedure to obtain the normalized strength function does not much affect the results. 
Indeed, we verified that the Gaussian smoothing instead of the Strutinsky smoothing leads to the essentially same 
results. 
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FIG. 7. (a)The partition function Xm(e) for m — 2 to 5 for the ISGQR corresponding to the normalized strength funciton 
plotted in Fig. 6 (a), (b) Its local scaling dimension _D m (e)for m = 2 to 5. The dashed curve represents D2(t) for the GOE. 



D. Origin of the Deviation from the GOE 



We shall investigate the origin of the observed energy scale C{ = 1.7MeV characteristic to the fine structures in the 
isoscalar quadrupole strength function. 



1. Doorway damping mechanism 



For this purpose we consider a damping mechanism by separating the shell model space to lplh and 2p2h subspaces, 
and we focus on the residual interactions within and between the subspaces. The giant resonance is a collective 
vibrational state whose dominant component is a coherent superposition of lplh configurations. The collective state 
can be described by the Tamm-Dancoff (TD) approximation, which corresponds to a truncated shell model where 
only lplh configurations are taken into account. The result of the TD approximation (TDA) is shown in Fig. 8(a), 
and is compared with the full shell model calculation including up to 2p2h's (Fig. 2(a)). It is seen that the peak 
position is well reproduced by the TDA. Note also that the isoscalar quadrupole strength in the TDA is spread over 
many eigenstates, which are correlated lplh states. This spreading of the quadrupole strength among lplh states is 
often called the Landau damping. Comparing Figs. 8(a) and 2(a), the Landau damping is found to give dominant 
contribution to the total FWHM of the strength function. 

If we take into account the 2p2h configurations, the interaction acting between lplh and 2p2h states (abbreviated 
as V12 ) takes part in. Since the 2p2h states have much more degrees of freedom, V\i causes spreading of the TD 
lplh states. Because of this spreading of lplh states, the quadrupole strength is also fragmented by Vi 2 on top of 
the Landau damping at the TDA level. To see the effect of Vi 2 , we performed a calculation where unperturbed 2p2h 
states and the interaction V\i is introduced. The resultant strength distribution, shown in Fig. 8(b), appears similar 
to the final result (Fig. 2(a)) as far as the gross profile of the strength distribution is concerned. 

Note also that there exists the residual interaction acting among 2p2h states, which we call V22 in the following. 
V22 causes configuration mixing among the 2p2h states. Comparison between Figs. 8(b) and 2(a) indicates effects of 
V22 on fine details of strength distribution, although the gross profile of strength distribution is not influenced very 
much. 

The above observations lead us to the following doorway damping picture. The quadrupole strength distribution is 
fragmented first in the lplh subspace. The lplh TD states then spread further over 2p2h states through V12. Here 
the lplh TD states may be considered as doorway states in the whole damping processes. The 2p2h states on the 
other hand may be regarded as states that play roles of background in the main damping processes since they do not 
dominate the total damping width, but they influence the fine structures through the coupling to the doorway states. 
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FIG. 8. (a) The strength function for the ISGQR obtained within the Tamm-Dancoff approximation and (b) the strength 
function for the ISGQR by omitting the coupling among the unperturbed 2p2h states. 

Assuming the above picture, we can consider energy scales that may be relevant to the fine structures of the 
quadrupole strength distribution. Concerning the lplh TD states which are considered as the doorway states, we 
have as energy scales 1) the spreading width (which we denote 712) of the TD lplh states caused by the interaction 
V12, in addition to 2) the average level spacing d\ p \h between the lplh states. Concerning the 2p2h states, 3) the 
level spacing d,2 P 2h between 2p2h states is to be noted. {d 2p 2h is almost identical to the level spacing d for the whole 
set of energy spectra since the number of 2p2h states is much larger than that of lplh states. ) Note also that 
mixing among 2p2h states caused by the interaction V22 is characterized by 3) the spreading width (722 named in 
the following) of 2p2h states. We have found previously that, in the case of a schematic model which incorporate 
doorway states coupled to background states, the spreading width of doorway states (corresponding to 712 in the 
present context) determines the energy scale where the local scaling dimension deviates from the generic GOE limit 
[14]. Note also that the deviation from the GOE is seen in the energy level fluctuation by the analysis of the A3 
statistics. The energy level fluctuation is sensitive to the residual interaction that causes mixing among unperturbed 
configurations. We found that the energy level fluctuation is governed by V22, but not by V12, as shown in Fig. 9. 
Thus we can consider another possibility that the energy scale of the deviation from the GOE limit may be related 
to the spreading width 722 of the 2p2h states. In the following, we shall focus on the effects of V12 and V22 which are 
connected to the two spreading widths 712 and 722 • 
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FIG. 9. The A3 statistics with rescaling the interaction strength V22 by a factor A22, for (a)A22 = 0.1 and (b)A22 = 2.0. The 
dashed curve represents the A3 for the GOE level fluctuation. We can find the large dependence of the A3 on the value of 
V22- We also verified that the A3 is insensitive to the change of the value of V12, namely the behavior of A3 at A12 =0.1 or 
A12 = 2.0 is almost same as that at A12 = 1.0, where A12 represents the rescaling factor of the interaction strength V12. 
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2. Spreading widths of lplh and 2p2h states 



Let us quantify 712 and 722- If we assume the Fermi golden rule, we obtain a simple estimate of 712 by 7f G = 

2 

2ir(lplh\Vi2\2p2h) /d,2 P 2h where (lplh\Vi2\2p2h) is a root mean square value of the matrix elements between lplh and 

2 

2p2h states. Similarly, we can make an estimation of the spreading width 722 as 7 2 = 2 G = 2tt {2p2h\V22\2p2h) /d,2 P 2h- 
The above estimate gives 7f G ~ 5.3MeV, and 7|2 G ~ 6.0MeV. Here we have inserted (lplh\Vi2\2p2h) = 98 keV, 
(2p2h\V22\2p2h) — 104 keV, d,2 P 2h — 11.3keV, which are evaluated for all the unperturbed lplh and 2p2h states in the 
energy interval of E = 20 — 30MeV. Note also di p ih ~ 0.5MeV for the lplh states in the same interval. Very similar 
values are obtained for (lplh\Vi2\2p2h) when we use the correlated lplh states (TD states) instead of unperturbed 
lplh's. 

It is also possible to evaluate 712 and 722 in a more direct way by using the strength functions of lplh and 2p2h 
states. Strength function of a lplh state \i, lplh) is given by 

Si,i p ih(e) = ^2 lplh\af 5{e - (E a - E hlplh )), (8) 

a 

where \a) and E a are the eigenstates of the full calculation and their energies, respectively, e denotes the relative 
energy from the energy centroid Ei_\ P \h of the lplh state considered. Evaluating an average of Si,ipi/i(e) over all the 
lplh states, the average strength function of lplh states is obtained, as shown in Fig. 10(a). Its FWHM gives an 
evaluation of the average spreading width 712 of lplh states. Similar method can be used to evaluate the spreading 
width 722 of the 2p2h states (cf. Fig. 10(b)). This direct method gives 7? 2 m = 1.5MeV and 7^ = 5.2MeV. 




E (MeV) E ( MeV > 

FIG. 10. (a) The strength function of the lplh state averaged over all lplh states within the energy interval E — 20 — 30MeV 
and (b) the strength function of the 2p2h states averaged over the same energy interval. 

It is noted that the values of the spreading width 712 are very different between the golden rule and the direct 
method. To investigate the difference, we evaluated root mean square values of interaction matrix elements of V12 for 
individual lplh states, and found that about one third of the lplh states have significantly large value ~ 120 — 200 keV 
as compared with the other lplh states which have small value ~ 30 — 70keV. Corresponding golden rule estimates 
of the spreading width are 7f 2 G = 8.3 - 22MeV and 7f 2 G = 0.5 - 2.8MeV, respectively. The FWHM of the average 
strength function of lplh states used in the direct method is dominated by the lplh states having small FWHM. 
Accordingly the direct method gives much smaller value (71™ = 1.5 MeV) than that of the golden rule (7f G = 5.3 
MeV). 



3. Dependence of LSD on V12 and V22 

In order to identify the origin of the observed energy scale ef = 1.7MeV, we shall analyze the dependence of the 
fine structures of the strength functions on V12 and V22 , which are intimately related to the spreading widths 712 and 
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722- For this analysis we rescale the interaction strength Vi 2 and V22 by factors A12 and A22 respectively, and look 
into the local scaling dimension to see how the fine structures change. 

It is seen from Fig. 11 that the fine structures probed by the local scaling dimension D m (e) strongly depend on A12. 
When the value of A12 is small, A12 = 0.1, the LSD exhibits very strong deviation from the GOE limit in almost all 
energy scales except at very small energy scale e/d <~ 1. As the value of A12 increases, the deviation becomes smaller 
and the LSD at small energy scales coincides with the GOE limit. The energy scale e* where the LSD start to deviate 
from the GOE limit moves to larger values as A12 increases. If A12 is increased more than the original value (e.g. 
A12 = 1.7), the LSD almost reaches the GOE limit over almost whole energy scales. The above dependence on A i2 
agrees very well with a similar behavior found in the simple doorway damping model [14]. (In that schematic model 
one can account for the spreading width 7d w of the doorway states analytically.) We found in the schematic model 
that the energy scale e* is proportional to 7d w as e* ~ 7dw/5 and the energy scale e** where the LSD decreases most 
steeply corresponds to 7d w - e* and e** moves with the interaction v between the doorway states and the background 
states since 7d w oc v 2 . This behavior is consistent with the A12 dependence in the present model. For quantitative 
evaluation, wc compare the spreading width 712 of the lplh TD states, which is considered here as the doorway 
states, and the energy scale ef where the deviation of the LSD from the GOE limit is maximum. Here we adopt the 
energy scale ef in place of e**, whose extraction is not easy in the present model due to fluctuation in the LSD. The 
energy scale ef = 1.7MeV found for the original interaction (i.e. A12 = 1) agrees with the spreading width 712 of the 
lplh states if we adopt the value 7^™ = 1.5MeV obtained by the direct evaluation. The agreement is a little worse 
if we refer to the value 7f 2 G = 5.3MeV derived from the golden rule estimate. We expect that the fine structures 
are dominated by the group of lplh states having smaller spreading width 712 ~ lMeV, while the remaining lplh 
states with 7 12 > lOMcV, which is larger than the total width of the giant resonance, do not contribute to the fine 
structure. The comparison with the spreading width jf™ by the direct evaluation would be more appropriate since 
7j 1 2 m reflects mainly the first group. (On the other hand, 7f 2 G based on the root mean square of interaction matrix 
elements puts emphasis on the other group. See the previous subsection.) The Ai2-dependence and the comparison 
with 712 indicate an approximate relation ef ~ 712, for which we estimate an ambiguity by a factor of two. 

Fig. 12 shows A22-dependence of the LSD. We immediately see that dependence on A 2 2 is much weaker than that 
on Ai 2 . For A22 < 0.2, the LSD exhibits deviation from the GOE limit for rather wide interval of energy scale, but 
the deviation itself is not very large. As A22 increases, the LSD approaches the GOE limit at small energy scales. (At 
A 2 2 = 0.5, for instance, the local scaling dimensions almost follow the GOE when e < lOd.) For A 22 > 0.5 the LSD is 
not sensitive to the change in A22 , and there remains the same deviation from the GOE limit even up to the maximal 
choice of A22 = 2.0. 

The insensitivity to A22 for A 2 2 > 0.5 may be understood as follows: We first note that the energy scale of the 
mixing among 2p2h states is given by the spreading width 7 22 of the 2p2h states, and that the mixing properties 
are expected to show up in the GOE behaviors for small energy scales satisfying e < 722 . In the case of the original 
strength A12 = A22 = 1, the spreading width 72™ =5.2MeV of the 2p2h states is larger than the spreading width 
Ti™ = L5MeV of the lplh TD states. If A22 < 0.5, opposite relation 722 < 712 is realized (n.b. 722 oc A| 2 and 
712 oc A^ 2 ), and hence the energy scale where the deviation from the GOE limit starts is determined by 722- In the 
case of A 2 2 > 0.5 (including the case of the original strength), on the other hand, we have 712 < 722, and the deviation 
from the GOE limit is governed by 712. In principle, we can consider a possibility that the LSD exhibits the two 
energy scales 712 and 722 separately. However the deviation related to 7 22 is small even when 7 22 < 712 (A 22 < 0.5), 
making it difficult to extract the energy scale associated with 7 22 . In the case of the original interaction strength, 
the spreading widths 7 2 ™ = 5.2MeV is close to the total width ^gr (~ 5MeV) of the giant resonance, and thus it is 
difficult to be probed. (Notice that we adopted a smoothing width 5MeV.) 

Summarizing the above two analyses, the characteristic energy scale of the fine structures is related to the spreading 
width 712 of the lplh TD states, which is considered as the doorway states in the giant resonance decay. 

4- Isovector strength function 

We have also analyzed the strength distribution for the isovector quadrupole operator. Figures 13(a) and (b) show 
the partition function and the local scaling dimension, respectively, for the isovector strength function in the energy 
interval E =20-30McV. The behavior of the LSD is similar to that of isoscalar strength function. Namely, the local 
scaling dimension almost follows the GOE limit at smaller energy scales whereas the deviation from the GOE limit is 
recognized at intermediate energy scales. We also verified that the change of the value of ef is quite sensitive to the 
change of A12 compared to the change of A22 as in the case of the ISGQR. Characteristic energy scale of the deviation 
reads ef ~ 70d ~ 0.9MeV. This value is smaller than that for the isoscalar strengths (ef <~ 1.7MeV) by about a factor 
of two. Both are consistent with the spreading width 7^™ = 1.5MeV. 
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(a) X, 2 =0.1 (b) ^2=0.4 




FIG. 11. The local scaling dimensions O m (e) (m = 2 to 5) for the ISGQR calculated by changing values of A12, which is 
chosen as (a) A12 = 0.1, (b)Ai2 = 0.4, (c)Ai2 = 0.7, (d)Ai2 = 1.0, (e)Ai2 = 1.4, and (f)Ai2 = 1.9. Here A12 represents the ratio 
between the adopted values of the Hamiltonian matrix elements between the TD lplh states and 2p2h states, and those of 
original values. The dashed curve represents .02(e) for the GOE. 
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(a) X 22 =0.1 (b) }. 22 =0.4 




e/d e/d 

FIG. 12. The local scaling dimensions D m (e) (m = 2 to 5) for the ISGQR calculated by changing values of A22, which is 
chosen as (a) A22 = 0.1, (b)A22 = 0.4, (c)A22 = 0.7, (d)A22 = 1.0, (e)A22 = 1.6, and (f)A22 = 1.9. Here A22 represents the ratio 
between the adopted values of the Hamiltonian off-diagonal matrix elements between 2p2h states and those of original values. 
The dashed curve represents f 2(e) for the GOE. 
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FIG. 13. (a)The partition function Xm( e ) (m = 2 to 5) for the IVGQR corresponding to the normalized strength funciton 
plotted in Fig. 6 (b). (b)Its local scaling dimension D m (e) for m = 2 to 5. The dashed curve rpresents D2(e) for the GOE. 

IV. CONCLUSIONS 

We have analyzed fine structures of the giant quadrupole resonances described by means of the microscopic shell 
model including up to 2p2h configurations. By applying the scaling analysis of strength fluctuations that utilizes 
the local scaling dimension, we extracted the energy scale e f = 1.7 MeV for the ISGQR and ef = 0.9 MeV for the 
IVGQR that characterize the fine structure. We discussed the origin of the characteristic energy scale in terms of the 
spreading widths of the lplh and 2p2h states. We found a clear correlation between the characteristic energy scale 
and the spreading width 712 of the lplh Tamm-Dancoff states, which play a role of doorway states in the present 
shell model description of the giant quadrupole resonances in 40 Ca. 

In the present paper, we intended to illustrate how the scaling analysis can be utilized to investigate the fine 
structures and their origin. Although the adopted model has some realistic features, we do not claim that the present 
shell model predicts all the quantitative features of the quadrupole strength distributions since the particle escaping 
width for instance is neglected in the present model. Also, other mechanisms such as the surface phonon coupling 
are not explicitly taken into account in the present model. Apparently one may need more realistic model if one 
intends to make a quantitative prediction or comparison with the experiment. Instead we would like to emphasize 
that the local scaling analysis provides a general tool for quantitative investigation of the fine structures, and the 
scaling analysis itself can be applied to any kinds of models which exhibit fine structures in the strength functions. 
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